Modulation of Superconducting Properties by Ferroelectric Polarization in Confined 

FE-S-FE Films 
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We show that the electric polarization at the interface with ultrathin superconducting (S) films 
sandwiched between ferroelectric (FE) layers allows achievement of substantially stronger modu- 
lation of inner carrier density and superconducting transition temperature as compared to FE-S 
bilayers typically used in superconducting FETs. We find that not only the larger penetration 
depths but also the pairing symmetry should be responsible for the fact that the electric field effect 
in high temperature superconductors is much stronger than in conventional systems. Discussing the 
advantages of multilayers, we propose a novel design concept for superconducting electric field-effect 
transistors based on ferroelectric films. 

PACS numbers: 74.80.-g,74.80.Dm,77.80.-e,79.60. Jv 



I. INTRODUCTION 



Electric field effect in superconductors continues to at- 
tract considerable attention in science and technology 1 . 
The external field can modulate the charge density 
and resistance, and control a reversible superconductor- 
insulator switching behavior which plays a key role in 
the development of superconducting field-effect transis- 
tors (SuFETs) 2 . Especially in the high-T c superconduct- 
ing cuprates, the field effect is expected to be strong, 
since their low carrier density leads to larger electric field 
penetration depths ~ 0.5 — 1 nm. Recently, complex fer- 
roelectric oxide materials with high dielectric constant 
like Ba^Sri^TiOa (BST) and Pb(Zr x Tii_ a: )03 (PZT)£ 
have been used to achieve substantial carrier modulation 
and shift of XL 5 of several K. 



Due to the small width of the accumulation layers at 
the contacts, the effect of the polarization-induced field 
is most significant in ultrathin superconducting films of 
a few nanometers thicknessi 3 ^. With the difficulties re- 
lated to the fabrication of ultrathin films and interfaces 
of good quality, theoretical modelling can be an addi- 
tional effective tool to study these systems. In the bilayer 
structures, the ferroelectric polarization at the interface 
attracts or repels the charge carriers in the supercon- 
ducting film, in close analogy to the effect of dopin g 2 : 3 : 5 . 
However, the field effect in FE-S multilayers where the 
polarization acts on both superconducting surfaces and 
may strongly modify the internal charge distribution in 
the entire ultrathin S-film, is a nontrivial and challeng- 
ing problem. Motivated by recent experiments^, in this 
work we show that the use of multilayers can drastically 
amplify the field effect, hence leading to much stronger 
increase of the superconducting transition temperature 
in comparison to S-FE bilayers typically exploited in 
SuFETs. 



II. THE MODEL 

In the approach considered here, we directly describe 
the internal charge redistribution in a ultrathin S-film 
caused by the ferroelectric polarization. Assuming the 
existence of an effectice attractive interaction, we focus 
on the question of how the superconducting state in the 
film can be influenced by this polarization. Specifically, 
we consider a hybrid system composed of periodically 
alternating FE- and S-layers with the ferroelectric po- 
larization P directed perpendicular to the interface, as 
shown in Fig. ^ The thickness of the FE subsystem is 
given by the number L\ of monolayers if = 1, . . . ,L\ 
in each ferroelectric layer. In the S-film, denotes the 
number of superconducting planes. We propose to de- 
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FIG. I: Scheme of the periodic FE-S multi-layer structure. 

scribe the structure by the following model: 
(i) The two possible orientations of ferroelectric dipoles 
(shown in Fig.QJ caused by the ion displacements on each 
lattice site (unit cell) / = 1, . . . , TV^ of the i/-th mono- 
layer in the FE-layer can be described by two values ±1/2 
of a pseudospin operator s(if,l). The spontaneous polar- 
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H MF 



H-£f}. The corresponding 



ization below the ferroelectric Curie temperature T C F can H % s s + H± 

then be modelled by an Ising-type Hamiltonian on the se t of the order parameters A is and P if is calculated self- 
cubic lattice^: Hp = — Jf s (*/j s (*/^')j where the 
summation ranges over all nearest neighboring sites; Jf 
denotes the internal dipolar interaction potential. The 
local polarization is given by the thermal average of Si f t i 
and is assumed to depend only on the monolayer index 



consistently by minimization of the system free energy. 
As a result, the proposed procedure describes the inner 
charge redistribution between the planes of an S-film in 
effective polarized medium on the boundaries. 



if- (si f ,i) = (s if ) = P if . 



(ii) In the superconducting film, we consider a BCS-like 
pairing in the planes i s — 1, . . . , L 2 : 
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The planes are coupled via the interplanar tunneling of 
the bosonlike pairs, 



H ± = /. Y^ih] h h, . , A . • h.C.) 



(2) 



the 



where the pair operators are bi a .k = Ci a ,—k,lCi s ,k,1: 
operator c\ k g creates an electron with spin a and im- 
pulse k in the i s -th plane, and Ck = A]£ (6 J k ). In @J, 

e k = —2t(cosk x + cos k y ) refers to the electron band en- 
ergy with lattice constant a = 1, and [i is the chemical 
potential. Assuming that some amount of charge N = 
2j k cr( n ia h o-) is either injected into the S-film from the 
electrodes or already exists due to chemical doping, we 
concentrate here on the effect of the carrier redistribution 
between the S-planes due to the contact polarization. For 
s-wave pairing, the superconducting pairing amplitude in 
each plane is: A ia = A k s = - J2 k ' V kk'( b is,k')' witn tlie 
effective pairing potential V?!, = — V la . 
(iii) At the interface we describe the screening of the sur- 
face charge due to the FE-polarization on the surfaces 
if = 1 and if = L\ by the electrons of the S-film. We 
consider here only the electrostatic interaction with the 
electron density of the boundary S-planes i s = 1 and 
i s = L2 given by the electron number operators rii a j , 
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Hi r , 



S1.1 ■ n L2 ,L. 



(3) 



where 7 is the charge-ferroelectric dipole interaction en- 
ergy. Note that, as the screening length in the cuprate 
superconductors is about 1 nm, the S-planes close to the 
boundary S-planes should also be taken into account in 
(O in the multilayered cuprates containing more than 
one Cu02 plane in the unit cell^ with a distance ~ 3.2 A 
between the planes. 

Furthermore, we focus here on the case where the 
FE-layers are far below the Curie temperature Tf . 
For this case, we neglect the polarization fluctua- 
tions and rewrite the energy Hp + Hi nt in the mean- 
field approximation, whereas the quasiparticle ener- 
gies for each given configuration |(1, ki, t), (l f — fei, I 
),..., (L 2 , fcL 2 j t)j {L2, —kL 2 i I)) are found by exact nu- 
merical diagonalization of the effective S-film energy 



n=0.1, j /t=_0._2 
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FIG. 2: Temperature evolution of polarization profiles in a 
heterostructure containing L2 = 2 S-monolayers for Jf /t = 1. 

The electron-dipole interaction 7 = eZuppj d 2 SF (e 
and Z are the electronic and the net FE-unit cell charges) 
depends on the distance dsF between the nearest FE- 
monolayers and S-planes. It depends also on the ionic 
displacement amplitude upp in the FE-unit cell. For 
instance, at dsF = 5 A, for a BST-layer with polar- 
ization P s ~ 26 /iC/cm 2 and ufe ~ 0.1 A, we obtain 
7 ~ 0.03 eV, whereas for ufe ~ 0.5 A(like in LiNbOa) 
we get 7 ~ 0.15 eV. A comparison with the case of a 
FE-S bilayer is easily provided by fixing the coupling 7 
at one of the contacts to zero. In our calculations, we 
take the pairing potential Vq = V la = 3.5t independent 
of i s . Considering the ideal surfaces, we should note that 
in SuFETs the surface roughness leads to the interface 
steps with a height h ~ 1 unit cell of SrTiC-3 as reported 
in Ref. 0- In the ultrathin S-films, such a step can act as 
a weak link and strongly affect the in-plane transport 10 . 
As long as h does not exceed the interplanar distance in 
the S-film, we can expect our results to not be signifi- 
cantly affected by the rough interfaces, but the problem 
of transport in this case needs more careful analysis. 



III. FERROELECTRIC POLARIZATION 

The selfconsistent solutions show stabilization of two 
different ferroelectric phases (+) and (— ) (discussed in 
details in Ref.|3) with decreasing temperature T, depend- 
ing on the electron concentration n = N / (N and 7. 
In phase (+) (Fig. [2J, the electric dipoles have the same 
orientation in each FE-layer (as shown in the inset of 
Fig-Bt a ))> an d the substantial deviations of the two possi- 
ble solutions for P (i^ 1 = — P^-i +1) due to the electron 
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screening @ appear close to the interfaces and decay ex- 
ponentially beyond the distance of about the ferroelectric 
correlation length 8 . As seen in Fig. [21 the polarization 
is also suppressed for smaller FE-layer thickness due to 
the depolarization field Ed = —JF^2((s(i't,l')) — Pbuik) 
acting on each dipole s(if, I), however the suppression ob- 
tained here is not as crucial as in Ref. and we still find 
a nonzero P for a layer with L\ = 6 monolayers. To get 
stronger suppression for a concrete FE-compound, one 
should consider more realistic long-range dipole-dipole 
interactions^ as well as strain effects. Nevertheless, we 
expect our main results related to the superconducting 
properties robust, since we consider the coupling 7 in the 
range 7/t < 1 (t ~ 0.1 eV for superconducting cuprates) 
which corresponds to the polarization P s < 25 /iC/cm 2 
obtained in Ref. ITU for FE-films of about 10 nm thickness. 
Also, taking into account the P-fluctuations, would sup- 
press the mean-field values of P considered in this work 
as an upper boundary for polarization profiles. With in- 
creasing coupling 7 (7/i > 1) and for higher electron den- 
sities n > 0.5, a transition into the ferroelectric phase (— ) 
occurs. The fluctuating polarization domains 7 appearing 
in the phase (— ), are a property of the Ising model in 
thin films. However, in a specific ferroelectric material 
these domains could be suppressed due to stronger depo- 
larization fields^. As the modulation of charge density 
in this state is much weaker than in phase (+), we focus 
here only on the polarization in the phase (+) which is of 
crucial importance for the superconducting properties. 



IV. REDISTRIBUTION OF CHARGE IN 
SUPERCONDUCTING FILM 

In the phase (+), the polarization P 1 > attracts the 
charge in the plane i s = 1 (Ani = ri\ — n > 0) while 
riL 2 decreases (Anj, 2 = nj, 2 — n < 0) as illustrated in 
Fig. |3f a). At low temperatures, the charge density redis- 
tribution due to P ^ 0, can be described by the difference 
An = rii — ul 2 ~ 3+ (-Pi + Pli) > 0. Consequently, the 
solution P 2 < leads to An < 0. Moreover, in con- 
trast to the bilayer system, the same direction of P in 
the second FE-layer, say, P > 0, repels the charge near 
the right contact and pushes it towards the accumula- 
tion region at the left contact, acting as an additional 
driving force for the charge supply into the accumulation 
plane i s = 1. This results in rii(multi) > n\(bi) and 
An(multi) > An(bi) which is clearly seen in Fig.^a). 

To understand the effect of the inter-planar charge re- 
distribution on the superconducting transition tempera- 
ture T c s , we consider the case L2 = 2. With the interpla- 
nar tunneling t± — > 0, the Hamiltonian Q)-© leads to 
two gap equations in the planes i s . For s-wave symmetry, 
different local Tf are found: 

T c s(1) oc 2iexp(-8i/V r + c5),T c 5(2) oc 2iexp(-8t/F - 5) ,(4) 

due to the factor S — g^fez^n (P1+P2) caused by the con- 
tact polarization. Hence, for low band filling (n < 1) and 
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FIG. 3: (a) Charge densities vs T in the S-planes i s = 1, 2 
in the multi- and bilayer systems containing Li = 10 FE- 
monolayers in the phase (+) with Jf/£ = 1, t±/t = 0.1. (b) 
The corresponding pairing amplitudes A ls are scaled by t. 



for P > 0, one finds that T C S(1) > T C S(2) and A 1 > A 2 , 
since the gaps A ls (T = 0) are comparable in energy to 

T c S(ia) . For Tc {2) < T < T C S(1) , this leads to a normal 
state in the plane i s — 2, whereas the plane i s = 1 is 
still superconducting. However, due to the intcrplanar 
coupling, even slightly above Tc^ 2 \ the superconduct- 
ing plane i s = 1 still induces a nonzero A 2 7^ for 



T < T, 



S(l) 



stabilizing the superconductivity in the en- 



tire S-film. This effect can be seen in Fig. Elb) from the 
numerical solution of the gap equations for t±/t — 0.1. 
For t± ^ 0, both gaps A 1 and A 2 vanish at a com- 
mon T = T c s which is higher than the local T C S ' 2 ' calcu- 
lated for t± = 0. The enhancement of the pairing on the 
planes without the pairing interaction solely due to the 
interplanar coupling is a well known property discussed 
in the literature^. In our case, however, the interpla- 
nar tunneling plays not only the role of an enhancement 
factor, but also provides a way for the accumulation of 
charge and stronger pairing in the accumulation S-planc 
caused by the polarization. In the multilayers, due to 
the stronger redistribution of the charge described by 
An, the last property leads to A 1 (multi) > A 1 (6i) and 
(multi) > T^Cbi) (Fig.[S{b)). For stronger coupling, 
the increase of T c s in the multilayer is also substantially 
stronger than that in the bilayer. The strongest increase 
of T c s caused by P can be achieved for low n; compared 
to the case 7 = (isolated the S-layer), we find a 25 — 35% 
increase of T c s with 7/i ~ 1 for n = 0.1, see Fig. 0J 



INTER-PLANAR COUPLING 



To estimate the effect of P for t± ^ 0, we rewrite the 
interplanar coupling in in the mean-field approxima- 
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tion and for Li = 2 introduce the order parameters as 

N±Al =-J2 V £'( b i,k')+Y, t * k 'M> ( 5 ) 
fe' fe' 

fe' fe' 

where the susceptibilities (6 Js ,fc) = g^fe- tanh (i/3-E£. s ), 

E^f are the eigenvalues of the mean- field superconducting 
Hamiltonian, and (3 — jr^p- The inter- layer tunneling 

1.24n 1 




FE-S coupling (y/t) 

FIG. 4: T c s vs coupling 7 in a system with L\ — 10 FE- 
monolayers for n — 0.1 and Jp/t = 1. 

t5* is introduced in a more generalized form so that 
the coherent tunneling (JSJ) corresponds to t± = tf* 5 kk > . 
As the field effect is expected to be stronger in high- 
T c superconductors, we consider further the two cases 
of s- and d-wave pairing with the ansatz A^ = A la rjk, 

v kk' = -V°VkVk' and<5i fc ' = t_L%%' in ©■ Here % = 1 
for the s-wave and rjk = cos(k x ) — cos(k y ) for the d-wave 
pairing symmetry. As for P > we have A 1 > A 2 in the 
case of isolated S-planes, we assume that the second term 
in the r.h.s. of the second equation (JSJ determining A 2 
is dominant and substitute it into the first equation for 
A 1 . The resulting quadratic equation has the following 
solution for small t±/V° -C 1 : 

T c s oc 2texp{-8t/V° + 8)), (6) 

which corresponds in fact to the upper local transition 
temperature given in J3J in the limit t± — > 0. However, 
in contrast to the expressions Q, here due to the inter- 
planar coupling we note three important features: (i) as 
analyzed in the literature, the pairing interaction is in- 
creased: V° = V°(l+(|^r) 2 ); (ii) the factor 6 = 

always leads to an increase of and does not depend 
on the direction of P, demonstrating that the enhance- 
ment of Xjf is essentially due to the S-planes with the 

dominant pairing (those with the higher local Tc in 
P]l): (iii) interestingly, the stronger inter-planar coupling 



makes the field-effect slightly weaker, since the contribu- 
tion of 5 in © decreases with the t± increase. The same 




coupling (y/t) 

FIG. 5: T c s versus 7 in P-S-P sandwiches and S-P bilayers 
for the cases of s- and d-wave symmetries calculated from (|SJ . 
Here n = 0.3, t±/t = 0.05, and the contact polarization Pi = 
Plx — tanh (|/3_E 9 ) fixed by an external field E g /t = 0.5. 

effect, i.e. the slightly weaker increase of T c s for larger 
t± is also seen in the numerical solution of iJSJ. The most 
striking fact, however, is that the coupling with the po- 
larization leads to a much stronger increase of T c s in the 
case of d-wave pairing when compared to s-wave, shown 
in Fig. [S] This result suggests that not only the larger 
screening length, but also the pairing symmetry should 
be responsible for the much stronger electric field effect 
in the high-temperature superconductors in comparison 
with the conventional compounds. Also, the difference 
between the T c s in the bilayers and sandwiches for the d- 
wave symmetry is much more pronounced which allows 
us to expect further progress in the field-effect experi- 
ments performed with FE-S-FE sandwiches. 



VI. SCHEME OF SANDWICH-BASED SUFET 

Having demonstrated the advantages of sandwiches, 
we propose a design of a field-effect transistor based 
on the confined FE-S-FE geometry. The power supply 
moves the opposite charge to the FE-gate electrodes. The 
operational principle of such an SuFET consists of two 
stages: (i) the superconducting state with the enhanced 
T c s can be reached in the regime where the system has 
the same direction of P in each FE-layer, as illustrated in 
Fig. H3 Note that for the fixed n, the antiparallel polar- 
ization in the FE-gates leads to lower charge densities in 
the accumulation region and consequently to the lower T c 
which is discussed in Ref. |F4J. (ii) To switch the SuFET 
into the insulating state, it is sufficient to destroy the 
accumulation layer at the first FE-S contact. This can 
be realized by switching the voltage V g — > — V g and tun- 
ing the polarization in the FE-layers to zero. Due to 
the simple geometry, even at the current stage of fab- 
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that in the proposed SuFETs also the effect of charge 
tunneling between the Cu02-planes could be tested in 
single- and especially in multilayered cuprates with inho- 
mogeneous carrier distribution 9 . 

In conclusions, we have shown that the FE-polarization 
in FE-S multilayers leads to much stronger modulation of 
inner carrier density and superconducting properties as 
compared to the bilayers. These advantages of multilay- 
ers can be used as the basis of novel design of proposed 
SUFETs. 



FIG. 6: Proposed scheme of SuFET based on confined FE-S- 
FE geometry. 
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